Weighted single step GWAS reveals genomic regions associated with economic traits in Murrah buffaloes.

The objective of the present study was to identify genomic regions influencing economic traits in Murrah buffaloes using weighted single step Genome Wide Association Analysis (WssGWAS). Data on 2000 animals, out of which 120 were genotyped using a double digest Restriction site Associated DNA (ddRAD) sequencing approach. The phenotypic data were collected from NDRI, India, on growth traits, viz., body weight at 6M (month), 12M, 18M and 24M, production traits like 305D (day) milk yield, lactation length (LL) and dry period (DP) and reproduction traits like age at first calving (AFC), calving interval (CI) and first service period (FSP). The biallelic genotypic data consisted of 49353 markers post-quality check. The heritability estimates were moderate to high, low to moderate, low for growth, production, reproduction traits, respectively. Important genomic regions explaining more than 0.5% of the total additive genetic variance explained by 30 adjacent SNPs were selected for further analysis of candidate genes. In this study, 105 genomic regions were associated with growth, 35 genomic regions with production and 42 window regions with reproduction traits. Different candidate genes were identified in these genomic regions, of which important are OSBPL8, NAP1L1 for growth, CNTNAP2 for production and ILDR2, TADA1 and POGK for reproduction traits.

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique.

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Effect of heat exposure on prostaglandin production and expression of COX-2, PGES, PGFS, ITGAV and LGALS15 mRNAs in endometrial epithelial cells of buffalo (Bubalus bubalis).

BACKGROUND: Early embryonic mortality is one of the major intriguing factors of reproductive failure that causes considerable challenge to the mammalian cell biologists. Heat stress is the major factor responsible for reduced fertility in farm animals. The present study aimed to investigate the influence of heat stress on prostaglandin production and the expression of key genes, including COX-2, PGES, PGFS, ITGAV and LGALS15, in buffalo endometrial epithelial cells. METHODS AND RESULTS: Buffalo genitalia containing ovaries with corpus luteum (CL) were collected immediately post-slaughter. The stages of the estrous cycle were determined based on macroscopic observations of the ovaries. Uterine lumens of the mid-luteal phase (days 6-10 of the estrous cycle) were washed and treated with trypsin to isolate epithelial cells, which were then cultured at control temperature (38.5 °C for 24 h) or exposed to elevated temperatures [38.5 °C for 6 h, 40.5 °C for 18 h; Heat Stressed (HS)]. The supernatant and endometrial epithelial cells were collected at various time points (0, 3, 6, 12, and 24 h) from both the control and treatment groups. Although heat stress (40.5 °C) significantly (P < 0.05) increased COX-2, PGES, and PGFS transcripts in epithelial cells but it did not affect the in vitro production of PGF2α and PGE2. The expression of ITGAV and LGALS15 mRNAs in endometrial epithelial cells remained unaltered under elevated temperature conditions. CONCLUSION: It can be concluded that elevated temperature did not directly modulate prostaglandin production but, it promoted the expression of COX-2, PGES and PGFS mRNA in buffalo endometrial epithelial cells.

Doppler ultrasonographic scan, gene expression and serum profile of immune, APPs and antioxidant markers in Egyptian buffalo-cows with clinical endometritis.

The objective of this study was to elaborate Doppler ultrasonographic scan, genetic resistance and serum profile of markers associated with endometritis susceptibility in Egyptian buffalo-cows. The enrolled animals were designed as; twenty five apparently healthy buffalo-cows considered as a control group and twenty five infected buffalo with endometritis. There were significant (p < 0.05) increased of cervical diameter, endometrium thickness, uterine horn diameter, TAMEAN, TAMAX and blood flow through middle uterine artery with significant decrease of PI and RI values in endometritis buffalo-cows. Gene expression levels were considerably higher in endometritis-affected buffaloes than in resistant ones for the genes A2M, ADAMTS20, KCNT2, MAP3K4, MAPK14, FKBP5, FCAMR, TLR2, IRAK3, CCl2, EPHA4, and iNOS. The RXFP1, NDUFS5, TGF-ß, SOD3, CAT, and GPX genes were expressed at substantially lower levels in endometritis-affected buffaloes. The PCR-DNA sequence verdicts of healthy and affected buffaloes revealed differences in the SNPs in the amplified DNA bases related to endometritis for the investigated genes. However, MAP3K4 elicited a monomorphic pattern. There was a significant decrease of red blood cells (RBCs) count, Hb and packed cell volume (PCV) with neutrophilia, lymphocytosis and monocytosis in endometritis group compared with healthy ones. The serum levels of Hp, SAA, Cp, IL-6, IL-10, TNF-α, NO and MDA were significantly (P˂0.05) increased, along with reduction of CAT, GPx, SOD and TAC in buffalo-cows with endometritis compared to healthy ones. The variability of Doppler ultrasonographic scan and studied genes alongside alterations in the serum profile of investigated markers could be a reference guide for limiting buffalo endometritis through selective breeding of natural resistant animals.

Oocyte competence and gene expression in parthenogenetic produced embryos from repeat breeder and normally fertile buffaloes (Bubalus bubalis) raised in sub-humid tropical climate.

The reproductive management of the buffalo species still faces several unresolved problems, which directly affect the productivity of the herd, one of them being the presence of repeat breeder females. Given this scenario, this study aimed to verify the developmental competence of oocytes obtained from repeat breeder females and submitted to parthenogenetic activation. In addition, embryo gene expression was compared to normally fertile females. Murrah buffaloes were divided into two groups: repeat breeder (RB, n = 8) and normally fertile or control (CR, n = 7). Cumulus-oocyte complexes (COCs) were aspirated by transvaginal ovum pick-up from estrus synchronized females. The COCs were submitted to IVM for 24 h, and subsequently, the oocytes were activated using ionomycin, followed by 6-DMAP. Afterwards, the presumptive parthenotes were cultured for six or seven days in a microenvironment of 5 % CO2, 5 % O2, and 90 % N2 at 38.5 °C. The expression of OCT4, GLUT1, BCL2 and TFAM genes from blastocysts was evaluated. The overall COCs recovery rate was 70.9 % (190/268). The maturation (57.8 vs 71.1), cleavage (45.2 vs 62.2) and blastocyst (30.1 vs 45.9) rates did not differ (P > 0.05) between RB and CR females, respectively. Similarly, no significant difference (P > 0.05) was observed for the expression of studied genes in both RB and CR females. In conclusion, oocytes obtained from RB were as developmentally competent as those collected from CR females, with similar energy metabolism and in vitro development capacity. Thus, the low fertility rate of repeat breeder buffaloes, when compared to normal cyclic females, must be due to subsequent events to the blastocyst stage.

Genome-wide 5'-C-phosphate-G-3' methylation patterns reveal the effect of heat stress on the altered semen quality in Bubalus bubalis.

Semen production and quality are closely correlated with different environmental factors in bovines, particularly for the buffalo (Bubalus bubalis) bulls reared under tropical and sub-tropical conditions. Factors including DNA methylation patterns, an intricate process in sperm cells, have an impact on the production of quality semen in buffalo bulls under abiotic stress conditions. The present study was conducted to identify DNA methylome signatures for semen quality in Murrah buffalo bulls, acclaimed as a major dairy breed globally, under summer heat stress. Based on semen quality parameters that significantly varied between the two groups over the seasons, the breeding bulls were classified into seasonally affected (SA = 6) and seasonally non-affected (SNA = 6) categories. DNA was isolated from purified sperm cells and sequenced using the RRBS (Reduced Representation Bisulfite Sequencing) technique for genome-wide methylome data generation. During the hot summer months, the physiological parameters such as scrotal surface temperature, rectal temperature, and respiration rate for both the SA and SNA bulls were significantly higher in the afternoon than in the morning. Whereas, the global CpG% of SA bulls was positively correlated with the afternoon's scrotal surface and rectal temperature. The RRBS results conveyed differentially methylated cytosines in the promoter region of the genes encoding the channels responsible for Ca2+ exchange, NPTN, Ca2+ activated chloride channels, ANO1, and a few structure-related units such as septins (SEPT4 and SEPT6), SPATA, etc. Additionally, the hypermethylated set of genes in SA was significantly enriched for pathways such as the FOXO signaling pathway and oocyte meiosis. The methylation patterns suggest promoter methylation in the genes regulating the sperm structure as well as surface transporters, which could contribute to the reduced semen quality in the Murrah buffalo bulls during the season-related heat stress.

Diversity assessment of a lesser known buffalo population from Central India and its comparative evaluation reveals presence of sufficient genetic variation and absence of selection.

Planned breeding and conservation strategies for a lesser-known population require an assessment of complete genetic diversity and population structure analysis in addition to its morphometric characteristics. In the present study, a comparative analysis of the genetic structure of a rare buffalo population, namely Chhattisgarhi, was extensively studied using a panel of FAO-recommended microsatellite markers along with well-established breeds namely Murrah, Nili-Ravi, Gojri, Kalahandi, and Nagpuri. Mode shift analysis indicated the absence of genetic bottleneck in the recent past. Assessment of genetic diversity indices across all loci indicated the presence of sufficient genetic variation within and between populations. Analysis of molecular variance between the six different buffalo populations attributed 19.05% of the variations to between-population differentiation. Cluster analyses using DAPC and Bayesian approach along with the phylogenetic tree based on UPGMA grouped six populations into three groups. The Chhattisgarhi population was revealed to be genetically closer to Nagpuri and Kalahandi populations. The study reveals the presence of sufficient genetic diversity within the Chhattisgarhi population and indicates the absence of a systematic selection program. We suggest improvement and conservation programs should be planned for this breed in the near future through short-term selection.

Exome-wide comparative analyses revealed differentiating genomic regions for performance traits in Indian native buffaloes.

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes.

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

A synonymous single nucleotide polymorphism (g.36417726C > A) in the Lama2 gene influencing fat deposition is associated with post-partum anestrus interval in Murrah buffalo.

Postpartum absence of estrus exhibition known as postpartum anestrus interval (PPAI) for more than 90 days after calving is a concerning issue for dairy buffalo farmers' economy. The PPAI duration is influenced by both management practices and animal genetics. Investigating genetic markers associated with PPAI is crucial for incorporating them into marker-assisted selection programs. Towards this goal, our study focused on exploring potential genetic markers from early postpartum adipose tissue gene networks. We successfully identified 24 Single Nucleotide Polymorphisms (SNPs) within 9 candidate genes. In our initial analysis involving 100 buffaloes, we detected a significant association (P = 0.02267) between a specific synonymous SNP within the Lama2 gene (g.36417726C > A) and PPAI. This finding was subsequently validated (P = 0.02937) in a larger cohort of 415 buffaloes, where the SNP explained 1.36 % of the genetic variance. Intriguingly, buffaloes with the CC genotype of this SNP exhibited a PPAI that was 12.71 ± 3.21 days longer compared to buffaloes with AA and CA genotypes. To gain insight into the functional relevance of this SNP, a computational analysis was performed which indicated that the C allele of the SNP (g.36417726C > A) increased the stability of LAMA2 mRNA compared to the A allele. This computational prediction was corroborated by observing a significant increase (P = 0.01798) in Lama2 gene expression (greater than 8-fold) and higher fat percentage (P < 0.05) in adipose tissue of CC genotypes (48.78 ± 1.87 %) compared to AA genotypes (33.59 ± 4.5 %). Furthermore, we noted a significant (P < 0.05) upregulation of C/ebpß, Pparγ, Fasn, C/ebpα, and Pnpla2 genes, along with the downregulation of Bmp2 and Ptch1 in CC genotypes as opposed to AA genotypes. This observation suggests the involvement of the Pparγ-mediated pathway in both adipogenesis and lipolysis within CC genotypes. In summary, our comprehensive analysis involving association and functional validation underscores the potential of the SNP (g.36417726C > A) within the Lama2 gene as a promising genetic marker against extended PPAI in Murrah buffalo.

Analysis of OCT4 and PGP9.5 gene expression in prenatal and postnatal buffalo (Bubalus bubalis) testes.

This study aimed to investigate and characterize the spermatogonial stem cells (SSCs) in buffaloes at different stages of development, including prenatal, neonatal, prepubertal, and adult testes. We sought a comprehensive understanding of these cells through a combination of histological, immunohistochemical, and ultrastructural analyses. Specifically, we examined changes in the expression of two potential SSC markers, OCT4 and PGP9.5, using immunohistochemistry. Additionally, we conducted a real-time quantitative polymerase chain reaction (RT-qPCR) to assess the relative gene expression of OCT4 and PGP9.5. The relative expression of the OCT4 gene was down-regulated in the adult testes compared to its expression during prepubertal and neonatal life. The relative expression of the PGP9.5 gene was up-regulated in the neonatal testes and down-regulated in the prepubertal and adult testes. The spermatogonia were round, oval-to-ellipsoidal cells lying over the basement membrane (BM) with a round-to-oval nucleus. Based on the immunoexpression of the putative SSC markers, OCT4 and PGP9.5, we concluded that the proportion of stem cells was highest during the neonatal stage, followed by the prepubertal and prenatal stages. This finding sheds light on the dynamics of spermatogonial stem cells in buffalo testes at different developmental stages, providing valuable insights into these cells' regulation and potential applications.

Identification of hub genes and their expression profiling for predicting buffalo (Bubalus bubalis) semen quality and fertility.

Sperm transcriptomics provide insights into subtle differences in sperm fertilization competence. For predicting the success of complex traits like male fertility, identification of hub genes involved in various sperm functions are essential. The bulls from the transcriptome profiled samples (n = 21), were grouped into good and poor progressive motility (PM), acrosome integrity (AI), functional membrane integrity (FMI) and fertility rate (FR) groups. The up-regulated genes identified in each group were 87, 470, 1715 and 36, respectively. Gene networks were constructed using up- and down-regulated genes from each group. The top clusters from the upregulated gene networks of the PM, AI, FMI and FR groups were involved in tyrosine kinase (FDR = 1.61E-11), apoptosis (FDR = 1.65E-8), translation (FDR = 2.2E-16) and ribosomal pathway (FDR = 1.98E-21), respectively. From the clusters, the hub genes were identified and validated in a fresh set of semen samples (n = 12) using RT-qPCR. Importantly, the genes (fold change) RPL36AL (14.99) in AI, EIF5A (54.32) in FMI, and RPLP0 (8.55) and RPS28 (13.42) in FR were significantly (p < 0.05) up-regulated. The study suggests that the expression levels of MAPK3 (PM), RPL36AL + RPS27A or RPL36AL + EXT2 (AI), RPL36AL or RPS27A (FMI) and RPS18 + RPS28 (FR) are potential markers for diagnosing the semen quality and fertility status of bulls which can be used for the breeding program.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes.

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Comparative RNA isolation methods from fresh ejaculated spermatozoa in Sahiwal cattle (Bos indicus) and Murrah buffalo (Bubalus bubalis) bulls for high quality and enhanced RNA yield.

Sperm mRNA transcriptional profiling can be used to evaluate the fertility of breeding bulls. The aim of the study was to compare the modified RNA isolation methods for higher RNA yield and quality from freshly ejaculated sperm of cattle and buffalo bulls. Ten fresh ejaculates from each Sahiwal (n = 10 bulls × 10 ejaculates) and Murrah bulls (n = 10 bulls x 10 ejaculates) were used for RNA isolation. From the recovered live sperm, total sperm RNA was isolated by conventional methods (TRIzol, Double TRIzol), membrane-based methods combined with TRIzol (RNeasy + TRIzol) with the addition of ß-mercaptoethanol (BME) and Kit (RNeasy mini) methods in fresh semen. Among different isolation methods; the membrane-based modified methods combined with TRIzol (RNeasy + TRIzol) with the addition of ß-mercaptoethanol (BME) resulted significantly (p < .05) higher total RNA quantity (300-340 ng/µL) and better purity in different concentrations of spermatozoa viz., 30-40 million, 70-80 million and 300-400 million sperm. The study concluded that the inclusion of BME to the combined membrane-based methods with somatic cell lysis buffer solution was best for constant increased yield and purity of RNA isolation from Sahiwal cattle and Murrah buffalo bull sperm.

Clonal spread of non-O157 Shiga toxigenic Escherichia coli O21:H25 in raw water buffalo milks.

AIMS: This study was conducted to investigate the presence of Shiga toxin-producing O157 and non-O157 E. coli in raw water buffalo milk, as well as to determine the virulence gene profiles, phylogroups, sequence types, and serotypes of the isolated strains. METHODS AND RESULTS: A total of 200 hand-milked raw water buffalo milk samples were collected from 200 different water buffaloes over a period of three months from 20 different farms. Isolation of STEC was performed using CHROMagar STEC. Presence of stx1, stx2, and eaeA genes were investigated by mPCR. Phylogroups and sequence types of E. coli strains were determined by Clermont phylotyping and MLST. Serotyping was performed using PCR or WGS. According to the results, two milk samples obtained from two different farms were found as STEC-positive. All Stx-positive E. coli isolates belonged to phylogenetic group A and were assigned to ST10. WGS results indicated that serotype of two isolates was O21:H25 and average nucleotide identity was detected at 99.99%. Thirteen additional registered E. coli O21:H25 assembled WGS data were obtained from EnteroBase and a phylogenetic tree was constructed. CONCLUSIONS: With this study, the presence of stx2 harboring E. coli O21:H25 in milk was identified for the first time. Although the identified serotype is considered a non-pathogen seropathotype, we conclude it could play an important role in the environmental circulation of Stx-phages and consequently contribute to the emergence of new STEC-related outbreaks.

Follistatin (FST) is expressed in buffalo (Bubalus bubalis) ovarian follicles and promotes oocyte maturation and early embryonic development.

Follistatin (FST), a member of the transforming growth factor-ß (TGF-ß) superfamily, has been identified as an inhibitor of follicle-stimulating hormone. Previous studies showed that it plays an important role in animal reproduction. Therefore, this study aims to investigate its effect on the maturation of buffalo oocytes in vitro, and the underlying mechanism of FST affecting oocyte maturation was also explored in buffalo cumulus cells. Results showed that FST was enriched in the ovary and expressed at different stages of buffalo ovarian follicles as well as during oocyte maturation and early embryo development. The FST expression level was up-regulated in MII buffalo oocytes compared with the GV stage (p < .05). To study the effects of FST on buffalo oocytes' maturation and early embryonic development, we added the pcD3.1 skeleton vector and PCD3.1-EGFP-FST vector into the maturation fluid of buffalo oocytes, respectively. It was demonstrated that FST promoted the in vitro maturation rate of buffalo oocytes and the blastocyst rate of embryos cultured in vitro (p < .05). By interfering with FST expression, we discovered that FST in cumulus cells plays a crucial role in oocyte maturation. Interference with the FST expression during the buffalo oocyte maturation did not affect the first polar body rate of buffalo oocyte (p > .05). In contrast, the location of mitochondria in oocytes was abnormal, and the cumulus expansion area was reduced (p < .05). After parthenogenetic activation, the cleavage and blastocyst rates of the FST-interfered group were reduced (p < .05). Furthermore, RT-qPCR was performed to investigate further the underlying mechanism by which FST enhances oocyte maturation. We found that overexpression of FST could up-regulate the expression level of apoptosis suppressor gene Bcl-2 and TGF-ß/SMAD pathway-related genes TGF-ß, SMAD2, and SMAD3 (p < .05). In contrast, the expression levels of SMAD4 and pro-apoptotic gene BAX were significantly decreased (p < .05). The FST gene could affect buffalo oocyte maturation by regulating the oocyte mitochondria integrity, the cumulus expansion, cumulus cell apoptosis, and the expression levels of TGF-ß/SMAD pathway-related genes.

Role of Different Members of the AGPAT Gene Family in Milk Fat Synthesis in Bubalus bubalis.

During triacylglycerol synthesis, the acylglycerol-3-phosphate acyltransferase (AGPAT) family catalyzes the conversion of lysophosphatidic acid to phosphatidic acid and the acylation of sn-2 fatty acids. However, the catalytic activity of different AGPAT members is different. Therefore, this study aimed to investigate the mechanism through which different AGPATs affect the efficiency of TAG synthesis and fatty acid composition. The conservation of amino acid sequences and protein domains of the AGPAT family was analyzed, and the functions of AGPAT1, AGPAT3, and AGPAT4 genes in buffalo mammary epithelial cells (BMECs) were studied using RNA interference and gene overexpression. Prediction of the protein tertiary structure of the AGPAT family demonstrated that four conservative motifs (motif1, motif2, motif3, and motif6) formed a hydrophobic pocket in AGPAT proteins, except AGPAT6. According to cytological studies, AGPAT1, AGPAT3, and AGPAT4 were found to promote the synthesis and fatty acid compositions of triacylglycerol, especially UFA compositions of triacylglycerol, by regulating ACSL1, FASN, GPAM, DGAT2, and PPARG gene expression. This study provides new insights into the role of different AGPAT gene family members involved in TAG synthesis, and a reference for improving the fatty acid composition of milk.

Electroporation-based CRISPR gene editing in adult buffalo fibroblast cells.

Electroporation is a widely used method for delivering CRISPR components into cells; however, it presents challenges when applied to difficult-to-transfect cells like adult buffalo fibroblasts. In this study, the ITGB2 gene (encoding the CD18 protein), plays vital for cellular adhesion and immune responses, was selected for editing experiments. To optimize electroporation conditions, we investigated parameters such as electric field strength, pulse duration, plasmid DNA amount, cuvette type, and cell type. The best transfection rates were obtained in a 4 mm gap cuvette with a single 20-millisecond pulse of 300 V using a 10 µg of all-in-one CRISPR plasmid for 106 cells in 100 µL of electroporation buffer. Increasing DNA quantity enhanced transfection rates but compromised cell viability. The 4 mm cuvette gap had high transfection rates than the 2 mm gap, and newborn cells exhibited higher transfection rates than adult cells. We achieved transfection rates of 10-12% with a cell viability of 25-30% for adult fibroblast cells. Subsequently, successfully edited the ITGB2 gene with a 30% editing efficiency, confirmed through various analysis methods, including T7E1 assay, TIDE and ICE analysis, and TA cloning. In conclusion, electroporation conditions reported here can edit buffalo gene(s) for various biotechnological research applications.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds.

BACKGROUND: Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated. RESULTS: Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds. CONCLUSION: The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Circulatory extracellular vesicle derived miR-195-5p promotes cellular apoptosis and suppresses cell proliferation in the buffalo endometrial primary cell culture.

In pregnant animals, communication between the mother and conceptus occurs via extracellular vesicles (EVs) that carry several biomolecules such as nucleic acids (miRNAs, mRNAs), proteins, and lipids. At the time of implantation, the endometrium undergoes several morphological and physiological changes, such as angiogenesis, apoptosis, and cell proliferation regulation at the implantation site, to attain a receptive state. This study was conducted to detect pregnancy-specific miRNAs derived from extracellular vesicles in the systemic circulation of Bubalus bubalis (water buffalo) and to assess their functional significance in the modulation of endometrial primary cells. The extracellular vesicles were isolated from the blood plasma using a precipitation-based method and further characterized by various methods such as Differential light scattering, Nanoparticle tracking assay, Western blot, and transmission electron microscopy. The relative expression of the selected extracellular vesicles associated miRNAs (EV-miRNA) at different intervals (days 15, 19, 25, and 30) post artificial insemination (AI) was analyzed using RT-qPCR, and expression of miR-195-5p was found to be significantly higher (P < 0.01) in pregnant animals on day 19 post AI (implantation window) as compared to day 15 post AI. The elevated expression might indicate the involvement of this miRNA in the maternal-conceptus cross-talk occurring during the implantation period. The KEGG pathway enrichment and Gene Ontology analyses of the miR-195-5p target genes revealed that these were mostly involved in the PI3-Akt, MAPK, cell cycle, ubiquitin-mediated proteolysis, and mTOR signaling pathways, which are related to the regulation of cell proliferation. Transfecting the in vitro cultured cells with miR-195-5p mimic significantly suppressed (P < 0.05) the expression of its target genes such as YWHAQ, CDC27, AKT-3, FGF-7, MAPK8, SGK1, VEGFA, CACAND1, CUL2, MKNK1, and CACAN2D1. Furthermore, the downregulation of the miR-195-5p target genes was positively correlated with a significant increase in the apoptotic rate and a decrease in the proliferation. In conclusion, the current findings provide vital information on the presence of EV miR-195-5p in maternal circulation during the implantation window indicating its important role in the modulation of buffalo endometrium epithelial cells via promoting cell death. Altogether, the milieu of miR-195-5p may serve as a novel and potential molecular factor facilitating the implantation of the early embryo during the establishment of pregnancy in buffaloes. Thus, miR-195-5p may be identified as a unique circulatory EV biomarker related to establishing pregnancy in buffaloes as early as day 19 post-AI.